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CO Core Candidates in the Gemini Molecular Cloud 


Yingjie Li 1 ’ 2,3 , Ye Xu 2 , Ji Yang 2 , Xin-Yu Du 2 ’ 3 , Deng-Rong Lu 2 , Fa-Cheng Li 2,3 

ABSTRACT 

We present observations of a 4 squared degree area toward the Gemini cloud obtained 
using J = 1-0 transitions of 12 CO, 13 CO and C 18 0. No C 18 0 emission was detected. 
This region is composed of 36 core candidates of 13 CO. These core candidates have a 
characteristic diameter of 0.25 pc, excitation temperatures of 7.9 K, line width of 0.54 
km s -1 and a mean mass of 1.4 Mg. They are likely to be starless core candidates, or 
transient structures, which probably disperse after ~10 6 yr. 

Subject headings: ISM: molecules - stars: formation - ISM: kinematics and dynamics 


1. Introduction 


Over the past several decades, relative to larger molecular gas surveys in Galactic plane, very 
few similar systemic studies were towards high g alactic latitudes. The Columbia surv ey conducted 
a large area survey up to a latitude b of ± 35° (jDame. Hartmann &; Thaddeu Il200ll h however its 
relative poor angular resolution (~ 8') is much larger than a typical core. 


The study of high-latitude molecular gas was active in the 1980s and 1990s, and had made 
some achievements: firstly, the high-latitud e molecular gas was researched by various of molecu les, 


such as CO and its isotopes, NH3, H2CO (IMag n ani. Blit z fe Wouter l oo t 198a : Turner 


1993 


a|0), 


etc.; secondly, in terms of its physical conditions, it is similar to the diffuse cloud in the galactic 
plane, and. in terms of it s ch emical abundance, it is similar to the cold dark cloud in the galactic 
plane (jTurner. Rickard &; Xu 1119891): and thirdly, t he high-latitude molecular gas was believed in 
the vicinity of sun (IMag nani . Blitz fe Mundy 19851): etc. 


We are carrying out a large-scale survey towards high-latitude molecular gas survey in the 
northern sky with a resolution of ~ 50". Our aims are to understand the structure, stability 
and physical conditions of molecular cores in the local area, between longitudes l from -10° to 
280° and b over ± 5°. For a starting, here we report our observations of the Gemini molecular 
cloud that is centered at l = 200° and b = 12°, and was also observed by the Columbia survey 
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(|Dame. Hartmann fc Thaddeu Il200lh . and no other literature made a further effort in this region 
in details. When it comes to “Gemini molecular cloud”, one usually associates it with “Gem OB 
1”. However, it has nothing to do with Gem OB 1, since its Vi sr is much less than Gem OB 1, and 
it is located at the east-north of Gem OB 1 with angular distance about 12°, and the distances to 
those two clouds are markedly different. With a resolution of 51 // , which is about a factor of 10 
higher than the Columbia survey, we are able to distinguish a molecular core with scale of about 
0.10 pc at a distance of 400 pc, which is estimated by combining several ways. See details in section 
3.2. 


2. Observations and Data Reduction 


We observed 12 CO (1-0), 13 CO (1-0) and C 18 0 (1-0) with the Purple Mountain Observatory 
Delingha (PMODLH) 13.7 m telescope from May 5 to June 21 and December 18 to December 
30, 2014. These three lines were simultaneously observed with the 9-beam superconducting array 
receiver (SSAR) working in the sideba nd separation mode and using the fast Fourier transform 
spectrometer (jShan. Yang &; Shi 2012!). 


Our observations were made in 16 cells of dimensions 30'x30', which covered an area of 4 
square degrees (248 pc 2 at distance of 400 pc). The cells were mapped using the_on-the-fly (OTF) 


obser vation mode with the standard chopper wheel method for calibration (IP enzias fe Burr us 


1 973 ). In this mode, the telescope beam scanned along lines of galactic longitude and galactic 
latitude at a constant rate of 50"/sec, and the receiver records spectra every 0.3 sec. Each cell was 
scanned in both the galactic longitude and the galactic latitude directions to reduce the fluctuation 
of noise perpendicular to the scanning direction. The typical system temperature ( T sys ) during 
observations was 250 K for 12 CO and ~160 K for 13 CO and C 18 0. Finally, we calibrated the 
antenna temperature (T 4 ) to the main beam temperature (T^) with a main beam efficiency (r^) 
of 46% for 12 CO and 51% for 13 CO and C 18 0. A summary of the observation parameters is listed 
in Table 1. 

Table 1 

Observation Parameters 


Line 

(J=l-0) 

(GHz) 

HPBW 

(") 

T 

- L sys 

(K) 

Vmb 

5v 

(km s 1 ) 

rms noise 
(K) 

12 CO 

115.271204 

49±2 

220-300 

46% 

0.16 

0.28 

13 CO 

110.201353 

51±2 

140-200 

51% 

0.17 

0.13 

C 18 0 

109.782183 

51±2 

140-200 

51% 

0.17 

0.19 
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3. Result 

3.1. General Distribution 

Figure 1 shows the distribution of 12 C0 emission in the Gemini molecular cloud. The strongest 
emission is located at (/, b ) = (200.26°, 11.57°) and (l, b) = (200.63°, 11.75°) with integrated 
intensities of 10.7 and 10.5 K km s _1 , respectively. Similarly, Figure 2 presents an integrated 
intensity map of 13 CO. The strongest emission is located at (l, b ) = (199.55°, 11.88°) and (/, b) = 
(199.95°, 11.78°) with integrated intensities of 1.1 and 1.0 K km s^ 1 , respectively. The positions of 
the emission peak of both lines are not exactly matched. Some emission in 12 CO has no counterpart 
in 13 CO , implying that the 12 CO is more widely distributed and extended than 13 CO. Unlike 
12 CO, 13 CO shows no continuous distribution, and only some condensations. Generally the gas 
distribution looks diffuse. Strong emission is only concentrated at several small locations. In the 
12 CO emission, clearly there are some filamentary structures that contain a few cores, while the 
13 CO emission is largely composed of separate, individual core. No C 18 0 emission was detected in 
a mean rms noise level of 0.19 K. We have averaged all of the spectra together to provide a single 
aggregate measure of the C 18 0 brightness over the cloud, which shows RMS of 0.01 K, and still 
unable to find any C 18 0 emissions. The reason may be that the emission of C 18 0 is too faint to 
be detected, and the upper limit for average Tj) of C 18 0 is 0.01 K. 

Figure 3 presents the longitude-velocity map of 12 CO. The longitude-velocity map of 12 CO 
suggests that V[ sr is confined to (-5, 5) km s _1 and mainly around ~ 0.0 km s _1 . A second velocity 
component (hereafter, velocity region 1) is centered at ~ 2.5 km s” 1 , ranging from l ~ 199.9° to 
l ~ 200.2°, and a third one (hereafter, velocity region 2) is centered at ~ -1.7 km s^ 1 in the range 
of l < 199.1°. There are large offsets from 0.0 km s^ 1 in the two velocity components. In addition, 
we find a velocity gradient of ~ 0.3 km s _1 pc -1 in the range of l ~ 200.0° to l ~ 200.7°. 

At first, to make clear of the velocity structure in those two velocity regions, the channel 
maps in Figure 4 illustrate the velocity structure of the molecular emission in part regions cor¬ 
responding to those two velocity regions. The left maps range from ( l,b ) = (199.89°, 11.70°) to 
(200.20°, 12.11°), while the right ones ranges from ( l,b ) = (198.79°, 12.10°) to (199.10°, 12.71°). 
Both sections correspond to velocity regions 1 & 2 in Figure 3. For ease of presentation, we refer 
to them as map 1 and map 2, respectively. 

Map 1 shows a montage of 12 CO emission distributions with velocities ranging from -2.5 to 5.5 
km s -1 in every 1 km s _1 , which indicates a drastic change of morphology: no cloud component 
appears in more than 2 successive channels. Similarly, map 2 shows the montage of 12 CO emission 
with velocities ranging from -3.5 to 2.5 km s _1 in every 1 km s _1 which presents a velocity dis¬ 
tribution of part of the filamentary components. This montage contains two distinct components 
of Vi sr that shift from -2.5 to 2.5 km s -1 . These complicated features in the velocity distributions 
likely indicate that the cloud consists of many components with velocities differing by ~ 5 km s _1 . 
Similar to map 1, map 2 also shows a big change of morphology; for instance, the bottom cloud 



-4- 



200.5 200.0 199.5 199.0 

Galactic Longitude 


Fig. 1.— Integrated intensity map of 12 CO , integrated between -7 and 7 km s 1 and scaled between 
0 and 11.0 K km s -1 . 

components appear in only 2 channels. The big change of morphology in the channel maps may 
be triggered by stochastic processes between clouds such as collisions and chaotic magnetic fields, 
rather than ordered motions such as rotation. Furthermore, unlike map 1, map 2 shows some 
filamentary structures, especially in the velocity range of -1.5 to -0.5 km s _1 and from 0.5 to 1.5 
km s -1 . In addition, 12 CO emissions seldom simultaneously appear at discontinuous channels in 
both maps. 

And then, we present the channel map of the entire region of the Gemini Molecular cloud in 
Figure 5. It shows that most 12 CO emissions are confined to (-2, 2) km s -1 , and indicates that the 
morphology in Figure 5 changes less than those two velocity regions shown in Figure 4. 


3.2. Distance Estimate 

In order to obtain the radius of a core candidate, denoted by R 13 for 13 CO, the mass (based 
on local thermodynamic equilibrium (LTE)), denoted by M, and the virial mass, denoted by My, 
the distance to the core candidate is indispensable. Kinematic distance is a widely used distance. 
Whereas the Vi sr of the cloud is around 0 km s -1 and cannot be effectively used to estimate it, 
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Fig. 2.— Integrated intensity map of 13 CO , integrated between -3 and 5 km s 1 and scaled between 
0 and 1.2 K km s _1 . The yellow circles show the positions of core candidates in 13 CO. 



Fig. 3.— Longitude-velocity map in 12 CO . 


therefore, we have tried other four methods to work it out. 
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Fig. 4.— Channel Maps of 12 CO in part regions. The two regions range from (l,b) = 
(199.89°, 11.70°) to (200.20°, 12.11°) in the left map, and from (l,b) = (198.79°, 12.10°) to 
(199.10°, 12.71°) in the right map. These sections corresponding to velocity region 1 & 2 in the 
longitude-velocity map. Each map was integrated over a velocity interval of 1 km s^ 1 , and the 
velocity ranges (in unit of km s -1 ) of each map are marked in the lower-left corner. The lowest 
contour is 3xRMS, and the contour interval is equal to (maxvalue of each map—minvalue of each 
map)/5. The core candidates of 13 CO fallen into the two regions are also indicated by red crosses. 
Note that the velocity range in both maps are the velocity range of the signal in each corresponding 
regions. 


Firstly, we try to find stars in the line of sight of the Gemini mol ecu lar cloud, a nd get a 
Tych o 2 star—TYC 1349-01421-1—with distance of 131 pc ( Ammons. Robinson. Strader et al. 
2006), which is located at (07:11:05.29, +17:00:45.9) in J2000 equatorial coordinate system. The 
Ak s = 0.174, and Ak = 0.95 Ak„ = 0.165, which indic ates that this star is s lightly reddened and 
may be also a foreground star as its Ak < 0.3 (see e.g., Lombardi. Alves fc Lada m3). 


Secondly, for a given cloud, the expectation value of the height above the plane is 0.798 a z , 
where a z is the Gaussian scale of the gas and is abou t 75 pc, and the distance to a cloud is 
0.798 a z esc 6, where b is the galactic latitude (see e.g.. IMaernani. Blitz fc Mundv 19851 ). Using this 
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Fig. 5.— Channel Maps of 12 CO in the whole regions. Each map is integrated over a velocity 
interval of 1 km s _1 , and the velocity ranges (in unit of km s -1 ) of each map are marked in the 
lower-left corner. 

method, we estimate that the distance is around 300 pc. However, this method is generally applied 
to clouds with b > 25°. Whereas, molecular gases are mainly in galactic plane. 

Thirdly, the facts that the Gemini molecular cloud is centered at l = 200° and the Vi sr is 
around 0 km s^ 1 indicate that this cloud could be nearby, for instance, like Taurns, Lam Ori 



















or Califirnia giant Molecular cloud whose distances are 140 pc, 400 pc and 450 pc, respectively 


(IJuvela. Lethtinen fe Mattila 1119971 : lLada. Lombardi fe Alves 1120091 : iLang. Masheder. Dame et al. 

2oooh . 


Lastly, we use the interstellar extinction distribution to estimate the distance t o Gemini molecu¬ 
lar cloud. We follow the approach applied by Marshall. Robin. Revie et al. (20061) and use th e stel¬ 


lar population synthesis model of the Galaxy constructed in Besangon (Robin, Reyle, Derriere et al. 


2003). Figure 6 shows the result. We find three 2MASS sources, J07084534+1610175, J07125749+1647429 
and J07104500+1755493, which are in line of sight of the 13 CO clumps of Gemini molecular cloud 
and have Ak s = 0.67(J — K s ) of 0.261, 0.526 and 0.532, respectively. And the distance estimated 
by 2MASS source J07084534+1610175 is 430T170 pc, with reliability of 0.68. However, this source 
has Ak = 0.95 Ak s = 0 .248 , which may indicate that this source may be a foreground star as its 
Ak < 0.3 (see e.g., Lombardi, Alves & Lada 2011), and therefore may slightly underestimate the 
distance. By using other two 2MASS sources which are definitely background stars, we find the 
distance reach to 2.6 kpc. While, both of those two 2MASS stars are about > 1 mag dimmer than 
J07084534+1610175, and have poorer Signal-to-Noise ratio with a factor of > 3.2. 

The Tycho 2 star—TYC 1349-01421-1 indicates a distance of 131 pc. The molecular gases 
mainly concentrate on galactic plane, and for molecular clouds of which b ~ 12°, the distance 
should be much less than 2.6 kpc; in addition, the second method shows a distance of ~300 
pc, which suggests that cloud is < 1 kpc away, and the third method also indicates that the 
Gemini molecular cloud is a nearby cloud. On the other hand, the forth method also has a large 
uncertainty. Furthermore, the CO emissions suggest an averaged column density of Hydrogen of 
0.7 xlO 21 cnr 2 (see details in section 3.4.4), which indicates that extinction in the J band_(Aj) is 
0.13 mag according to the relation of JVh 2 = 5.57 x 10 21 Aj cm~ 2 ma g -1 (IVuong et al. 2003 ), and 
the Ak s ~ 0.1 by using extinction curve of Aj/Ak s = 0.282/0.114 ( Cardelli . Clay ton. &; Mathi s 
19891 ). The Ak s ~ 0.1 suggests a distance of ~ 120-170 pc, and also supports the cloud is far less 


than 1 kpc away, but the uncertainties of this relation are the scale coefficient and extinction curve, 
which are difficult to estimate errors here. And the Ak 3 ~ 0.1 derived from the averaged column 
density of hydrogen may give an evidence that source J07084534+1610175 is most likely to be a 
background star. Combining those factors, we therefore adopt the distance of 400 pc, and denote 
it by d re f. Because the distance estimated from the first method and averaged column density 
of Hydrogen is less than 200 pc, this distance may be overestimated. The units of any quantities 
depended on distance would multiply by an additional factor— d/d re f = d/(400 pc), where d is in 
unit of pc, which is the genuine distance to the Gemini molecular cloud. 


3.3. Core Identification 

The core candidate catalogue produced in this paper is created from the 13 CO data cube. We 
define cores as overdensities in a molecular cloud that may or may not contain cores from which 
single or multiple stars are born. These are also known as star-forming or starless cores, respectively 
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Fig. 6.— Result of distance estimation by interstellar extinction distribution, left: J — K s .vs. 
distance, each asterisk represents a bin from our method. Right: corresponding J — I\ s histogram 
of the model and the observations. The black line represents the 2MASS observations and the red 
line is the result of our method. 


( Parsons. Thompson. Clark et ah 2012; Williams. Blitz & McKee 2000) 


In order to decompose the Gemini molecular cloud into cores, we looked to automated core 
detection programs to deal with the obtained 3-dimensional FITS cube file of 13 CO. In this pa- 


CLUMPFIND (Williams, de Geus & Blitz 

1994), which has been widely used in the literature 

(Moore, Bretherton, Fuiivoshi et al. 

2007; 

Buckle, Curtis, Roberts et al. 

2010 

). Due to the in- 

trinsic biases that are present in automated core detection routines ( 

Schneider &; Brooks 

2004; 


Smith. Clark &; Bonnell 2008), we chose the core candidates carefully, while advising caution when 
performing blind cross-comparison of the properties of the core candidates measured here with 
those from other core detection algorithms. The CLUMPFIND algorithm first contours the data 
and searches for peaks in order to locate the core candidates, and then follows them down to lower 
intensities. To obtain as much of the emission as possible without contamination from noise, we set 
the parameters TLOW=2xRMS and DELTAT=2xRMS in 13 CO, where TLOW determines the 
lowest level contour of a core candidate, and DELTAT represents the ga p betwe en contour levels th at 


det ermines the lowest level at which to resolve merged core candidates (Williams, de Geus & Blitz 


1994). The parameters of each core candidate, including the position, velocity, size in the galactic 


1 http://st arlink.jach.hawaii.edu/starlink 
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longitude and galactic latitude directions, and one-dimensional velocity dispersions, are directly ob¬ 
tained using this process. Moreover, further steps were taken to ensure the detection of real cores. 
First, we excluded core candidates that had voxels touching the edge of the map area in order to 
avoid core candidates where the full extent of the emission may not have been recovered. Secondly, 
we excluded core candidates whose sizes were smaller than the beam resolution. Thirdly, we ex¬ 
cluded core candidates that had ratios of two main-axes greater than 4. Finally, the morphology 
of the core candidates were checked by eye within the three-dimensional galactic longitude-galactic 
latitude-velocity space to identify core candidates with meaningful structures. Using these criteria, 
we obtain 36 core candidates (plotted in Figure 2 and catalogued in Table 2) in the 13 CO maps. 

Table 2 

Properties of the 36 core candidates of 13 CO in the Gemini molecular cloud 


Core name 

Vlsr 

(km s -1 ) 

Ai>13 
(km s _1 ) 

Ris 

{{d/d def ) pc) 

Tex 

(K) 

T13 

at h 

(xl0 21 cm -2 ) 

M 

((d/d def ) 2 M 0 ) 

M v 

«d/d def ) M 0 ) 

G198.81+12.29 

-2.32 

0.40 

0.05 

7.7 

0.77 

1.4 

0.3 

0.9 

G198.87+12.22 

1.66 

0.41 

0.11 

7.7 

0.20 

0.4 

0.3 

2.4 

G198.88+12.64 

-1.83 

0.69 

0.13 

7.3 

0.25 

0.7 

0.7 

7.3 

G199.04+12.13 

1.00 

0.33 

0.07 

7.4 

0.59 

0.8 

0.3 

0.9 

G199.35+11.72 

1.66 

0.29 

0.05 

7.4 

0.68 

0.8 

0.1 

0.5 

G199.51+11.81 

0.17 

0.77 

0.18 

10.7 

0.09 

0.5 

1.0 

13.3 

G199.55+11.87 

1.49 

0.62 

0.13 

9.9 

0.30 

1.3 

1.3 

5.8 

G199.57+11.57 

1.49 

0.63 

0.11 

6.9 

0.31 

0.7 

0.7 

5.5 

G199.72+11.80 

0.17 

0.68 

0.09 

8.2 

0.25 

0.9 

0.7 

5.8 

G199.74+11.87 

-0.17 

0.50 

0.15 

8.0 

0.23 

0.5 

0.7 

4.4 

G199.81+12.12 

-0.83 

0.45 

0.07 

6.9 

0.39 

0.7 

0.3 

2.0 

G199.84+11.09 

0.17 

0.50 

0.18 

9.1 

0.26 

0.8 

1.7 

5.6 

G199.95+11.77 

3.98 

0.42 

0.09 

8.7 

0.36 

0.8 

0.7 

2.2 

G199.97+11.89 

3.49 

0.63 

0.07 

9.5 

0.27 

1.1 

0.3 

3.3 

G199.98+12.07 

3.32 

0.13 

0.07 

8.4 

0.50 

0.3 

0.1 

0.2 

G200.01+12.07 

2.49 

0.48 

0.07 

8.4 

0.30 

0.7 

0.3 

2.2 

G200.03+11.90 

2.99 

0.25 

0.09 

9.5 

0.34 

0.6 

0.3 

0.7 

G200.14+11.90 

0.50 

1.60 

0.25 

8.2 

0.16 

1.3 

5.3 

80.2 

G200.15+12.22 

-0.33 

0.85 

0.11 

6.8 

0.35 

1.1 

1.0 

9.8 

G200.18+12.37 

-0.83 

0.81 

0.18 

7.7 

0.25 

0.9 

2.3 

15.5 

G200.20+11.57 

0.00 

0.47 

0.13 

8.3 

0.27 

0.6 

0.7 

3.3 

G200.22+12.16 

-1.66 

0.52 

0.11 

7.3 

0.28 

0.6 

0.7 

4.0 

G200.25+12.25 

-0.33 

0.62 

0.09 

7.1 

0.22 

0.5 

0.3 

4.5 

G200.26+11.57 

-0.66 

0.70 

0.11 

7.6 

0.26 

0.8 

0.7 

6.9 

G200.26+12.11 

-1.16 

0.57 

0.18 

7.2 

0.24 

0.5 

1.0 

7.1 

G200.29+11.53 

0.66 

0.63 

0.05 

7.6 

0.28 

0.8 

0.3 

3.1 

G200.29+11.04 

-1.00 

0.35 

0.05 

8.4 

0.39 

0.7 

0.1 

0.9 

G200.31+12.21 

-1.16 

0.37 

0.11 

6.7 

0.36 

0.5 

0.3 

2.0 

G200.37+11.49 

-1.33 

0.43 

0.07 

7.5 

0.37 

0.7 

0.3 

1.5 

G200.42+11.77 

-0.17 

0.58 

0.20 

8.2 

0.23 

0.7 

2.0 

8.7 

G200.46+11.47 

1.00 

0.47 

0.36 

8.3 

0.31 

0.7 

6.9 

10.4 

G200.47+11.98 

0.17 

0.49 

0.36 

7.4 

0.32 

0.7 

5.6 

10.7 

G200.58+11.77 

-0.33 

0.17 

0.07 

7.6 

0.54 

0.4 

0.1 

0.2 

G200.61+11.76 

-1.00 

0.41 

0.09 

7.6 

0.44 

0.8 

0.3 

1.8 

G200.63+11.82 

0.50 

0.60 

0.09 

7.6 

0.33 

0.9 

0.3 

4.2 

G200.69+11.77 

-0.17 

0.54 

0.09 

7.0 

0.31 

0.6 

0.3 

3.1 

mean 

0.32 

0.54 

0.12 

7.9 

0.33 

0.7 

1.1 

6.7 
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Properties of the 36 core candidates of 13 CO are summarized in Table 2. The columns represent, 
respectively, the catalogue core (candidate) name, Vi sr , the line width (Awi 3 , be., the velocity 
component along the line of sight, which is fitted by Gaussian profile), its radius (-R 13 ), its excitation 
temperature ( T ex ), the optical depth (T 13 ), the column density of hydrogen molecular (-/Vh 2 ), its 
LTE mass (M) and its virial mass (My). Core candidate names also convey their location, for 
instance: G199.04+12.13 represents the galactic longitude 199.04° and galactic latitude 12.13°. 
The last row in Table 2 is the mean value of each physical quantity. 


3.4. Physical Properties 

3-4-. 1. Radius 

The core candidate radii are derived from the geometric mean of the core candidate sizes in 
two directions, and we have made the size of the beam deconvolved from the radius measurement. 
Figure 7 shows the distributions of the core candidate radii in 13 CO. The mean core candidate 
radius in 13 CO is 0.12 ( d/d re f ) pc, and the median core candidate radius of 13 CO is 0.10 (d/d re f) 
pc, which is slightly less than the mean value. 
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Fig. 7.— Distribution of core candidate radii in 13 CO (in unit of (d/d re f) pc). The mean and 
median values (in unit of ( d/d re f ) pc) are marked on the upper-right corner of this plot. 
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3-4-2. Excitation Temperature 


To obtain the excitation temperature of each core candidate, three assumptions are required. 
First, the core candidates are in local thermodynamic equilibrium (LTE); secondly, the emission 
line used in this determination is optically thick; and_thirdly, the emission fills the beam. Assum¬ 


ing T ex ( 13 CO) = T ex ( 12 CO) = T ex ( see e.g.. Keto fc Mvers 1986), we calculated the excitation 
temperatures of the 12 CO emission using the following: 


T P , 


hv 




/ fcr*( 12 co) i 

\ hv exp (hv/kT^g) — 1 



-i 


( 1 ) 


where T ex is in unit of K, hv is the energy of a single photon emitted during the transition of 
12 CO(J=1-0), k is the Boltzmann constant, T R ( 12 CO) is the peak main beam temperature in unit 
of K (for 12 CO emission), and T^ g is 2.7 K, i.e., the temperature of cosmic microwave background 
radiation. 

Excitation temperatures of the 36 core candidates catalogued in Table 2 range from 6.7 to 10.7 
K, with a mean temperature of 7.9 K. We note that the excitation temperatures may be under¬ 
estimated when we consider the effects of beam dilution and CO self-absorption on the obtained 
observations. 


3-4-3. Opacity 

The opacity of 13 CO is given by 


T* r { 13 CO) 

5.29 (1/ (exp (5.29/T ex ) - 1) - 0.164) 

where T ex and T^( 13 CO) (in unit of K) are given in Section 3.4.2. This formula indicates that 
uncertainties in the excitation temperature will directly affect the opacity. The opacities of the 36 
core candidates are catalogued in Table 2. We have calculated the ratio of T^( 12 CO)/TJ l ( 13 CO) in 
the position of peak emission of those 13 CO core candidates, and this ratios range from 3.1 to 7.2, 
with mean ratio of 4.2, and we then find that the 12 CO is indeed optically thick with t ±2 > 1, and 
the 12 CO opacity of 58.4% core candidates are larger than 5. 




(see e.g., Kawamura, Onishi, Yonekura et al. 1998): 


3-4-4- Column Density 


The column density of 13 CO 1V( 13 C0) is given by (see e.g.. K awamura . Onishi . Yon ekura et al. 


1998) 
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iV( 13 CO) = 2.42 x 10 14 x 


”M3^ f ! ;r A 7 ; 13 

1 - exp (-5.29 /T ex ) ’ 


(3) 


where the units of iV( 13 CO) are cm“ 2 , T ex is the excitation temperature of the J=l-0 transition of 
13 C0 in unit of K. The line width AU 13 is in unit of km s . The value of T ex , T 13 and AU 13 are 
catalogued in Table 2. 


Column densities in ter ms of H? as opposed to 13 C0 are determin ed using the abundance 

ratio between H 2 and 13 C0 ([Wilson Hl999l : lMagnani. Blitz fc Mundv Ill985l ). and the value adopted 

here is 1.2xl0 6 . Ah 2 of the 36 core candidates catalogued in Table 2 range from 3.4xl0 20 to 
, —2 


1.4xl0 21 cm -2 with a mean value of 7.5xl0 2U cm " 2 (see table 2). Figure 8 shows the distribution 
of Ah 2 . Systematic errors in the column densities derived here may arise from both their opacity 
and their abundance ratio. 



Column Density 


Fig. 8 .— Distribution of column densities (in unit of xlO 21 cm 2 ) of molecular hydrogen. The 
mean and median values (in unit of xlO 21 cm -2 ) are marked in the upper-right corner. 

Surface densities in unit of g cm -2 are calculated by multiplying the column densities by the 
mass of molecular hydrogen. Within the observed Gemini molecular hydrogen, the surface densities 
range from 1.5xlCP 3 to 6.4xl0 ~ 3 g cm" 2 , with mean value of 3.4xl0 ~ 3 g cm" 2 . 
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3-4-5. Mass 

Once the density and size of a core candidate have been determined, the LTE mass (M, in 
unit of (d/d re f) 2 M 0 ) of the core candidate reads 


M = ANi3qq 


[h 2 ] 

FcoJ wnH ’ 


(4) 


where A = 7r(i?i 3 ) 2 is the area of the core candidate in unit of cm 2 , where R\s (in unit of (d/d re f) 
pc) is the 13 C0 core candidate radius which is given in Table 2 and section 3.4.1. Ni 3 0 o is 
the column density of 13 C0 in unit of cm -2 , as given in section 3.4.4., and the abundance ratio 
is again assumed as [H 2 ]/[ 13 CO]~ 1.2 x 10 6 . /in 2 is the mean molecular weight of gas per H 2 
mol ecular (mh-. = 2.72 that includes hyd rogen, helium and the isotopologues of carbon monoxide; 


e-g-, 


Buckle . Curtis . Roberts et al. 2010 ) and mu is the mass of a single hydrogen atom. 


The core candidate masses estimated here, which are based on the assumption of LTE, range 
from 0.1 to 6.9 (d/d re f) 2 M 0 with mean mass of 1.1 (d/d re f) 2 M 0 (see Table 2). The mass error 
comes from the distance error and its abundance ratio. Because some CO emission is absorbed by 
grains, these estimates can be considered as lower limits, and an estimate of the intervening dust 
mass needed to be investigate to determine how much it will affect our observations. The total 
mass of these 36 core candidates is 38.3 (d/d re f) 2 M 0 . 

To obtain the virial mass My, the radius and line width are required, and a core density profile 
of p oc r~ 2 is assumed. The virial mass My reads: 


My = 126i?i 3 (Au) 2 , 


(5) 


where My is the virial m ass in unit of (d/ddef) M 0 , Rva is the radius in ( d/d re f ) pc, and Av is the 
line width in km s _1 (MacLaren, Richardson & Wolfendale 19881). The uncertainties in the virial 
masses come from the measurement errors of line widths and i?i 3 . 


The virial masses derived for 13 CO are about 7.0 d re f/d times larger than the LTE masses, 
and have a range of 0.2 to 80.2 ( d/dd e f) M 0 . The a mean virial mass is 6.7 ( d/ddef) M 0 (see Table 
2 ). 


3-4-6. Vi sr and Line Width 

The core candidate line widths increase with increasing LTE masses, however the excitation 
temperatures do not appear to have this tendency. This implies that thermal motion is not the 
exclusive broadening mechanism due to (JThermai ~ ( T/A) 3 / 2 km s” 1 oc T 1 / 2 in a coarse sense, 
where (JThermai is the thermal line width (line width broadened via pure thermal motion), T is the 
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kinetic temperature in unit of K, which > T ex , and A is the atomic mass number. 

To determine the dominant broadening mechanisms, we calculate thermal line widths of tracer 
species in unit of km s -1 using aThermal = y/kT/m/ 1000, where k is the Boltzmann constant, m 
is the mean molecular mass in unit of kg and T is the kinetic temperature, which > T ex , in unit of 
K. The non-thermal line width o n on-Thermai and thermal line widths of gas as a whole <JThermai,g 
read: 

& Non—Thermal = \J °\d ~ a Thermal = \] 13 / “ ° Thermal’ (®) 


1 / kT 

^Thermal,g — i nnn 1 / 

1000 y /iWH 


(7) 


where am is the one-dimension al v elocity dispersion in unit of km s^ 1 , and AU 13 is the li ne width 
in unit of km s -1 (see e.g., Myers 1 983 ); and k is the Boltzmann constant, /r = 2.4, mn is the 
mass of a single hydrogen atom in unit of kg and T is the kinetic temperature, which > T ex . in 


unit of K. aThermal ranges from 0.04 to 0.06 km s -1 , (JThermal,g ranges from 0.15 to 0.19 km s _1 and 
the distribution of a^on-Thermai is plotted in the top histogram of Figure 9. Because T > T ex , 
the aThermal and aThermal,g is underestimated (i.e., a lower limit) and a^on-thermai is therefore 
overestimated (i.e., a upper limit). 


The mean ratio of a^on-Thermai to aThermal,g is 1-35 (see the bottom histogram in Figure 
9), which implies that the non-thermal broadening mechanism plays a dominant role in the core 
candidates; or to be more details, 19% (i.e. 7/36) of those are subsonic core candidates, and the rest 
are supersonic core candidates, in which the non-thermal broadening mechanism plays a dominant 
role. 


The core candidates of high Vi sr (greater than 2 km s^ 1 ) and low V) sr (lower than -1.5 km s -1 ) 
in 13 CO belong to two different regions (i.e., map 1 and map 2 respectively) which are plotted in 
the channel maps seen in Figure 4. 

The big change of morphology in the channel map in Figure 4, imply that the velocity gradients 
mentioned in section 3.1 (i.e., 0.3 km s _1 pc -1 ) and the large dispersion of Vi sr of the 13 CO core 
candidates (i.e., 6.3 km s _1 ) in the Gemini molecular cloud may be more strongly affected by 
stochastic processes, such as collisions and chaotic magnetic fields, rather than ordered motions 
such as rotation. 


4. Discussion 

In this section, we will discuss the statistical properties and star-formation processes occurring 
in the Gemini molecular cloud. There is no sharp demarcation between these two aspects. In fact, 
statistical properties can be regarded as tools in research of star-formation. 
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Fig. 9.— Top: histogram of non-thermal line widths. The mean value and median value are 
marked in the right-upper corner of the histogram. Bottom: histogram of non-thermal (blue) line 
width and thermal (red) line width of gas as a whole. We numbered the 13 CO core candidates in 
ascending order of galactic longitude, and marked them on the x-label for a brief view. 

4.1. Core Candidate Mass Function 

The core mass function (CMF) describes the relative frequency of cores with differing masses. 
CMF is commonly fitted with a basic power law dN/dM oc M - “, where N is the core candi- 
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dates in each bin, M is the mass and a is the corresponding power index. The shape of the 
CMF has been seen to resemble the stellar initial mass function (IMF) that descri bes the rel ative 


frequen cy o f stars w ith differing masses over a large range of environments (e.g., Salpeter 1955.; 


Simpson, Nutter &; Ward-Thompson 1120081 ). Under standing the relationsh i p between the CMF an d 


the IMF can help constrain star formation models (Bate fc Bonne ll 2005|; 


Reid & Wifson 2006) 


To obtain the parameter index of the CMF, we fixed the bin widths and counted the number 
of core candidates per bin. The value of abscissa is the mean value in each bin, and two kinds of 
CMF derived from 13 CO are given by: 


dN 

dM 


oc M 


-ai 


( 8 ) 


dN 

dMy 


oc My a \ 


(9) 




Fig. 10.— CMF (where mass includes M in the units of (d/ddef ) 2 M 0 and My in unit of ( d/dd e f) 
Mq). The black star-like points are logarithmic values of the number of core candidates per unit 
mass against logM (left) and log My (right). a\ and a .2 are 1.83±0.12 and 1.63±0.12, respectively, 
cci and «2 correspond to power the index of CMF of LTE mass and viral mass as derived from 
13 CO, respectively. 

where M and My are the LTE masses (in unit of (d/d re f ) 2 Mq) and virial masses (in unit of 
(d/d re f) Mq) derived from 13 CO, respectively. N is the number of core candidates in each bin, 
while ot\ and «2 are the corresponding power indices. The power law correlations of CMF are: 
dN/dM = (0.57 ± 0.04)M- L83±ao7 (left), and dN/dMy = (0.99 ± 0.11)MR L63±0 12 (right), see 
Figure 10. Note, the significance of the coefficients of 0.57 ±0.04 and 0.99 ±0.11 is very limited for 
the uncertainty of distance. Our power index (1.83 for LTE mass and 1.63 for virial mass derived 

















- 18 - 


from 13 CO , see F igu re 1 0) is bet ween the IMF_of 2.35 ( 
presented by Kramer, Stutzki, Rohrig et ah (1998) and 


e.g. , Salpeter 1955 ) and CMF of J^.6 -1.8 


Buckle, Curtis, Roberts et ah (2010) 


The similarity between the CMF and IMF power indices can be simply be explained by a 
constant star formation efficiency that is unrelated to the mass and self-simila r clo ud structure, and 
based on a scenario of one-to-one tra nsformation from cores to stars (ILada. Muench. Rathborne et al. 
2008). In addition, a simulation by Swift & Williams (12 00 8) suggested that the obtained IMF is 
very similar to the input CMF even when different fragmentation modes are considered. However, 
the full physical meaning of the relationship between the CMF and the IMF remains unclear due 
to a number of complications, including: (a) completeness limitations, (b) time-scales, (c) ph ysical 


size scales, (d) distances, etc. (Curtis & Richer 2010; Reid, Wadsley, Petitclerc et ah 2010) 


4.2. Star Formation 


Figure 11 shows the relationship between the virial mass My and LTE mass M. The virial 
parameter ay , defined as My/M , describes the competition of internal supporting energy against 
the gravitational energy. Virial parameter is inversely proportional to distance, and its typical value 
is 7.0 (d re f/d), and all virial parameters for the 36 core candidates in 13 CO are larger than 1.9 
( ddef/d )• The massive core candidates tend to have lower virial parameters. Furthermore, a linear 
fitting between LTE mass ( M in unit of (d/d re f) 2 M 0 ) and virial mass (My in unit of ( d/d re f) M 0 ) 
may be a good indicator of the typical virial parameter: the first-order coefficient is 4.84 (ddef/d ), 
as shown in Figure 11, which is slightly less than the typical virial parameter value of 7.0 (d re f/d). 
However, the correlation coefficient of this linear correlation is just 0.60. 

The mass relationship can also be fitted with a power-law of My oc M~< (where M in unit of 
(d/d re f) 2 M 0 and My in unit of ( d/d re f) M 0 ), and the result is My = (0.76±0.15)M a9 ' ±o ' 15 , c.c. = 
0.85, where c.c. is the correlation coefficient, and the significance of the coefficient of 0.76 ± 0.15 is 
very limited for the uncertainty of distance. 


The power index value of 0.97 that we obtained is lar ger t han the value of 0.67 in Orion B (which 
are gravitationally bound) reported by Ikeda fc Kitamura (20091) and 0.61 in Planck cold clumps 
(which are gravitationally bound) presented by Liu. Wu &; Zhangl ( 2012), and is significantly highe r 


than the index of pressure-confined clumps (ay oc M 2 / 3 ) given by Bertoldi Sz Mckee (1992). 


While the virial parameters are larger than 1.9, which indicates that the core candidates may be 
unbound. 


Figure 11 also shows that the virial mass changes more dramatically in the lower end of the 
LTE mass than in the upper end of the LTE mass. Additionally, the ratio between the virial and 
LTE masses decreases as LTE mass increases. 

However, the discussion of whether these core candidates are pressure confined or not may be 
not necessary if such structures do not need to be long lived. We then calculate a time-scale (t;) 
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Fig. 11.— Virial mass My (in unit of ( d/d re f ) M 0 ) .vs. LTE mass M (in unit of ( d/d re f ) 2 M 0 ) 
relation of the core candidates. The blue slanted line: the slope is 1.0, denotes the minimum value 
of virial parameter. The red slanted line: the slop is 4.84, and the correlation coefficient is 0.60, 
presents the linear correlation among virial mass and LTE mass. The green curve: shows the power- 
law correlation among virial mass and LTE mass, and the result is My = (0.76 ± 0.15)Af°- 97±0 ' 15 , 
c.c. = 0.85, where c.c. is the correlation coefficient. 


used by lLarson (119811 ') : 


ti ~ 4.6 x 10 6 


Rl3 


d 


2 y/2 In 2cr/ d re f 


yr, 


( 10 ) 


where, Ri 3 is in unit of pc, and a t = y 7 ( aNon-thermai) 2 + {(^Thermal, g ) 2 is the overall velocity 
dispersion of the mean particle in the cloud, where a t , a n on-thermal and (JThermai,g are in unit of 
km s -1 . The typical value of r/ is about 1.7 x 10 6 yr. And thus, the Gemini molecular cloud may 
be transient structure. 


Cores t hat do not h ave a k nown infrar ed associations a re starless cores or pre-protostellar cores 


(Benson & Mvers 

1989; 

War cl- T h o m n s o n. Scott. Hills et al. 

1994) 


1994). Those cores represent a stage 


earlier than protostellar. We have checked the IRAS point source catalogue in this region and found 
that no 13 CO core candidates are associated with IRAS point sources. We also have checked the 
WISE 22 /im data, and found that only one core candidate (i.e. Core candidate G200.15+12.22) is 
associated with WISE source J071257.49+164743.5, which has instrumental profile-fit photometry 
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magnitude of 9.0 mag in the band of 22 /jm. However, the columns of w4sigmpro = “null” and 
w4chi = 0.9, which indicates that this source is not measurable in the band of 22 /mi, where the first 
two letters of “w4” in each column mean band 4 with wavelength of 22 /mi. We then conclude that 
all 36 13 CO core candidates may be starless core candidates. Nevertheless, to determine whether 
these core candidates are in the process of forming stars, or already formed stars, or even other 
situations, further study is necessary. 


5. Summary 

We presented PMODLH mapping observations for an area of 4 deg 2 toward the Gemini molec¬ 
ular cloud centered at l = 200° and b = 12° in 12 CO, 13 CO and C 18 0 lines. The main results are 
summarized as follows: 

1. We identified 36 core candidates in 13 CO. Derived from 13 CO, we yield typical radii, column 
densities and LTE masses of 0.12 ( d/d re f ) pc, 7.5xl0 20 cm~ 2 and 1.1 (d/ddef ) 2 Mg respectively. 
We also found the mean excitation temperature of 7.9 K as derived from 12 CO. The total LTE 
mass derived from these 36 core candidates in 13 CO is 38.3 ( d/d re f) 2 Mg. 

2. Non-thermal broadening mechanism plays a dominant role in 81% (29/36) of those core 
candidates. 

3. Velocity gradients (i.e., 0.3 km s^ 1 pc -1 ) and large dispersions in Vi sr of the 13 CO core 
candidates (i.e., 6.3 km s _1 ) in the cloud may be more largely affected by stochastic processes, such 
as collisions, and chaotic magnetic fields rather than ordered motions such as rotation. 

4. Two kinds of CMF are present in this paper, namely the CMF of LTE masses, and the 
virial masses derived from 13 CO. Their power indices are a\ = 1.83 and ct 2 = 1-63, respectively. ot\ 
and «2 are all lower than the power index (2.35) of the IMF. Moreover, the CMF of LTE mass is 
steeper than the CMF of the virial mass derived from 13 CO. 

5. 13 CO core candidates in the Gemini molecular cloud are more likely unbound, and massive 
core candidates have lower virial parameters. All 36 core candidates are potential starless core 
candidates, and they also may be transient structure. 
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We also acknowledge Xiaolong-Wang, Chong Li and ShaoBo Zhang for their valuable help. This 
work is supported by the National Natural Science Foundation of China (grant Nos. 11133008 and 
11233007), the Strategic Priority Research Program of the Chinese Academy of Sciences (grant No. 
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